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'Cluster' extensions of the dynamical mean field method to include longer range correlations are 
discussed. It is argued that the clusters arising in these methods are naturally interpreted not as 
actual subunits of a physical lattice but as algorithms for computing coefficients in an orthogonal 
function expansion of the momentum dependence of the electronic self-energy. The difficulties with 
causality which have been found to plague cluster dynamical mean field methods are shown to be 
related to the 'ringing' phenomenon familiar from Fourier analysis. The analogy is used to motivate 
proposals for simple filtering methods to circumvent them. The formalism is tested by comparison 
to low order perturbative calculations and self consistent solutions. 
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I. INTRODUCTION 

Over the last decade the 'dynamical mean field' 
method^ has emerged from earlier investigations- of the 
infinite dimensional limit as a very useful tool for theoret- 
ical investigation of correlated electron systems. It pro- 
vides (within a certain approximation) a nonperturbative 
means of obtaining the electron self-energy and spectral 
function, and allows treatment of inelastic and thermal 
effects on the same footing as ground state energetics. It 
has revealed new insights into the physics of the Mott 
transition^ of 'colossal' magnetoresistance manganites?i 
of plutonium^ and many other systems, and may be com- 
bined with LMTO band theory (for reviews see RefsiH 
to treat correlation effects in realistic models. 

As originally formulated^ the dynamical mean field 
method is a local approximation. The problem which one 
actually solves is a single-site model self-consistently em- 
bedded in a medium. The results are an exact representa- 
tion of the physics of lattice models only in a limit of infi- 
nite coordination number. 1-2 While the method captures 
local physics very well, the treatment of intersite correla- 
tions is an important open issue. Attempts to formulate a 
controlled expansion about the infinite coordination limit 
have not led to useful and tractable expressions. An al- 
ternative approach focuses on the self-consistent embed- 
ding of a larger cluster in a medium^Si^iiSiiiiiSiiii^ and is 
closely related to attempts to extend the 'single-site CPA' 
approximation to more than one siteii^*i& Most of the 
approaches published to date involve choosing a specific 
cluster (subset of sites of the actual lattice of interest) 
as well as a specific embedding (geometry of connection 
of cluster sites to sites of the medium). An alternative 
('E-DMFT') method involves using hybrid boson-fermion 
methods to treat two-particle intersite correlationsi2ii2iiS 
and has had some success. 

While these approaches have led to a number of inter- 
esting results (for recent examples see e.g. Refs. ITolEI 
l2ll) . one worries that the choice of specific geometry both 
of the cluster and its embedding in the medium may bias 



the physics. Further, the real-space-cluster-based meth- 
ods sometimes lead to self energies which are non-causal 
for some momenta and frequencies. This is regarded as a 
grave defect of the method and its cure is an important 
open problem. 

In this paper we put forward an alternative point of 
view on the question of extending dynamical mean field 
theory to include intersite correlations. Our motivation 
is as follows. The original local dynamical mean field 
method, although often described in terms of a site self 
consistently embedded in a medium, may also be de- 
scribed as the 'momentum-independent self-energy ap- 
proximation'. It amounts to 

(i) replacing the general position-dependent self- 
energy S(i?, uj) by the on-site value S(0, u), and 

(ii) providing a prescription for computing E(0,w) 
from the solution of a single-site quantum impurity 
model. 

Considering the problem more generally, one may write 
£(p, u) as an expansion in orthogonal functions and trun- 
cate the expansion at a finite order, thereby replac- 
ing a general function by a small number of frequency- 
dependent coefficients, which may be determined from 
the solution of a several site fictive impurity model: "fic- 
tive" because although the several-site impurity model 
can be regarded as a cluster self-consistently embedded 
in a medium, this cluster need not be a sub-cluster of 
the physical lattice considered. It should be viewed as 
merely a device for computing the self-energy functions 
of interest. The general idea of taking an abstract view 
of the 'impurity' or cluster in the dynamical mean field 
method is not new. It is mentioned in Ref. and 
has recently been elegantly exploited by Kotliar and co- 
workers to link band theory and the dynamical mean field 
method. 6 It is also the basis of an approach to the prac- 
tical solution of the single-site DMFT equations put for- 
ward by Potthoffi^ Here we apply the idea to the study 
of spatial correlations. 

The idea that the cluster model is merely an algorithm 
for computing coefficients in a low-order orthogonal func- 
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tion expansion of a physical self energy gives insight into 
the 'causality violations' which sometimes occur in dy- 
namical mean field schemes. In the context of cluster ex- 
tensions of dynamical mean field theory, the term 'causal- 
ity violation' means that for some momenta and frequen- 
cies, the imaginary part of the approximate lattice self- 
energy computed from the cluster approximation has the 
'wrong' sign. We shall argue that this is nothing but 
the 'ringing' phenomenon familiar from Fourier analysis: 
if one approximates a non-negative but sharply peaked 
function by a few low order terms in an orthogonal func- 
tion expansion, the approximant will change sign. We 
shall use the analogy to propose several simple cures. 
We n ote that the 'DCA' formalism used for example in 
Refs. I8I9H21I corresponds to a choice of orthogonal func- 
tions which are everywhere non-negative but which have 
no common support. In this sense, 'ringing' is absent, 
so the method is plainly causal, but discontinuities in 
momentum space occur. 

The outline of this paper is as follows. In section II 
we outline the formalism we shall use. In section III we 
present examples of the use of several-site impurity mod- 
els to calculate the self-energy of the two dimensional 
Hubbard model. Section IV discusses 'causality viola- 
tions' and possible cures. Section V is a conclusion. 



II. FORMALISM 

We study a model of electrons moving on a d- 
dimensional lattice with short ranged hopping ampli- 
tudes (described by a quadratic Hamiltonian Hq) and 
interactions (described by a Hamiltonian Hint)- Physi- 
cal properties of the model may be derived from the gen- 
eral 'Luttinger-Ward' expression for the action, which we 
may write in terms of the exact Green function G of the 
model as 



S = Trln(-G) 



kel 



fGl 



(1) 



with f&sfeei [G] the sum (with appropriate symmetry fac- 
tors) of all vacuum to vacuum 'skeleton' diagrams drawn 
with full Green functions (G) and no self-energy inser- 
tions. The electron self energy X may be obtained via 



£ = S$ skel /6G. 



(2) 



The functional S is defined for any G once the inter- 
actions (Hi n t) are specified. The correct G for a given 
'band structure' (-Ho) is determined from the equation 
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(3) 



It is convenient for our purposes to eliminate G in favor 
of X via a Legendre transformation (this transformation 
is also basic to the work of Potthof 22 ) defining 



Note that it follows from Eqs. lll'l 1[) that 

5$skel 



G = — 



ST, 



(5) 



and that in this representation the theory is fixed by de- 
manding that the G[S] obtained from Eq. JHJl is identical 
to the G obtained via 



G[S1 = Gn'-E 



(6) 



i.e. by minimizing with respect to X the functional 

fi[E] = -Tr In (-Gq 1 + S) + $ sfeei [£]. (7) 

The original momcntum-indcpendent-self-energy 
(single-site) dynamical mean field approximation may 
be formulated from Eqs. 1 1 1( > B as follows (essentially 
this derivation is given in Ref. 1): define $i oc as the 
approximation to the exact ^ s hel [S(p, u))] obtained by 
replacing the exact momentum dependent self-energy 
by a local approximant S; oc (lj) = j(dp)H(p, u>) which 
depends only on frequency. Consistency demands that 
this be equivalent to replacing $ s kei [G] in Eq. by the 
analogous quantity defined with the local Green func- 
tion. Because S/ oc (w) depends only on frequency, the 
quantity which follows from Eq. JSJ is the local Green 
function G/ oc = J(dp)G(p,uj) and Eq. © becomes the 

relation G/ oc = J(dp) [Gq 1 (p.w) - £; oc (w)] . 

The crucial observation which makes the single-site 
dynamical mean field approximation useful is that be- 
cause $i oc is a functional only of a function of frequency, 
it may be defined nonperturbatively as the solution of 
a single-site (quantum-impurity) model which is speci- 
fied by a frequency-dependent Weiss field and by interac- 
tion terms related to the local interactions of the original 
model. The Weiss field is fixed by demanding that Eq. 
© is satisfied when the left hand side of this equation is 
the Green function calculated from the impurity model 
and the right hand side is the local Green function cal- 
culated from the lattice Hamiltonian, using the impurity 
model self-energy. 

A generalization is now evident. Consider a set of func- 
tions {4n, ipi} such that 



so that 



E(p,w)=5^(p)£<(w) 



(8) 



(9) 
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Tr (EG) ; 



(4) 



with Ej(u;) = J (dp) tpi(p)'E(p, us). Inserting Eq. JjJ) into 
Eq. (gj yields 

$ sfcei [{SJ] = $ sfce; - Tr ( S < G *) ( 10 ) 

i 

with G. t = / (dp) &(p)G(p,a;). 



3 



Now, construct an approximant for the self-energy as 
a sum of a small number of terms in the expansion given 
in Eq. © 



J approx 



few) 



E 



i(p)Si(w). (11) 



cluster to the medium must be addressediZii&iS while 
in the 'DCA' the Laue-function arguments advanced by 
Aryanpour et al. 1A require a momentum-independent in- 
teraction. A separate paper will analyze the issue from 
the present point of view^ 



Define & a pprox[^approx] as the functional obtained from 
<&skei by using the approximate self-energy instead of 
the exact X. Eq (JSJ) implies that this construction is 
equivalent to approximating & s kei by the set of dia- 
grams drawn using only the G conjugate [in the sense 
defined below Eq. (|10|) ] to the retained Ej. We see 
that ^approx is a functional of n + 1 frequency depen- 
dent fields. It therefore corresponds to the solution to 
some (n + l)-site fictive impurity model involving n + 1 
Weiss fields, and interactions derived (as discussed be- 
low) from the original model. The Weiss fields are fixed 
by the requirement that the impurity-model Green func- 



tions Gi 



c/#Sj are given by appropriate in- 



tegrals over the lattice Green function: 



G t (u) = /(dp)^( P )[Go 1 (p,w) 



From this point of view we may interpret the impurity 
model simply as a mathematical means for calculating, 
nonperturbatively, an approximation to the self-energy. 

The original dynamical mean field method corresponds 
to retaining only the i — term in the self energy. In 
the "DCA" approach of Jarrell and co- workers^ the func- 
tions 4>i are obtaining by tiling the Brillouin zone into re- 
gions TZi and setting 4>i (p) = 1 if p is contained in TZi and 
<pi (p) = otherwise. These functions are clearly nonneg- 
ative everywhere, and are orthogonal because they have 
no common support: at any p exactly one (pi is non- 
zero. This choice of functions leads approximants with 
discontinuities in momentum space. Other choices are 
discussed below. 

To completely specify the impurity model we must de- 
termine the interaction terms. Specifying the interaction 
terms is a subtle issue in general. One approach is to 
observe that the skeleton functional, and therefore its 
approximant Q a pprox is defined for any hopping Hamil- 
tonian H n . We may therefore consider the special case 
of no hopping (so H is simply the energies of what- 
ever on-site levels are considered). For models (such as 
the Hubbard model) in which the interaction is local, 
both the Green function and self-energy are diagonal in 
the site representation and are easy to compute. Com- 
parison of the exact and approximate expressions then 
shows (as was already demonstrated for the single-site 
DMFT in Ref. 0) that the interaction terms in the im- 
purity model are simply the original interaction terms of 
the lattice model. However, for longer ranged interac- 
tions, the situation may become more complicated. In- 
deed the difficulty with longer range interactions appears 
to be common to all schemes. For example, in real space 
cluster schemes the issue of interactions connecting the 



III. EXAMPLE APPROXIMANTS 

A. General considerations 

A cluster extension of dynamical mean field theory 
involves finding an impurity model to represent the 
frequency-dependent expansion coefficients in Eq. J^J. 
Such a model must involve n fields which have an or- 
thogonality property, so that we may unambiguously de- 
termine n independent Green functions and self energies. 
One convenient choice is to introduce an m > n compo- 
nent spinor of fermions tp and write an action of the form 

S = S {2) + S mt (13) 
with (we have suppressed the spin indices here) 



b * M * 



.i=0,. . . ,m 



(14) 



where bi are frequency dependent Weiss fields and the 
m x to matrices Mi satisfy 

Tr[M i M j ]=mS ij . (15) 

A general impurity model Green function is given by 

G lmp =goMo+g 1 M 1 +g 2 M 2 + ... (16) 

with coefficients gi(u>) given by 



1 5\nZ, 
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imp 



777 Sbi 

The orthogonality relations imply the self energies 



(17) 



m 



Mi 



Q 1 

imp 



(18) 



leading to self consistency equations of the form Eq. 112(1 . 
The self consistency condition in the local (no-hopping) 
limit implies Mq is the unit matrix. 



B. Harmonic expansion— second order 

Here we write, for a hypercubic lattice in d dimensions 
Vapproxip^) = E„M + ^e^ a E a (w) (19) 
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with a = ±x, ±y... so that naively there are 2d + 1 mean 
field equations: 



Go = / (dp)G p (u), 



G a 



(20) 
(21) 



so we should write an impurity model action involving 
2d + 1 fields. However, we may argue that in a cubic 
lattice, all of the components S a are equal, so that the 
physics may be expressed via an impurity model depend- 
ing on two fields, So and S a . There are two self consis- 
tency conditions. The impurity model is then specified 
by a partition function arising from a functional in- 
tegral over the action 



S (2) = (V4 V: 



b Q b% 
h b 



S 



(2) 



(22) 



The mean field equations may be written in sym- 
metrized form as 



2Sb 
2Sh 



(dp)G p (tu), 
(dp) l P G p (uj) 



where 



G p {w) 



1 



e p — E — 2d7 P Si 



(23) 
(24) 

(25) 



One may consider including longer ranged 
correlations-for example, in d > 1 the second neighbor 
correlation, writing 

^ ap prox(p,^) = Sp (a;) + y~V p ' a I]i(tj) 

a 

+ e t{p - a+p - b) Y, 2 {uj). (26) 

We must then seek a multi-site impurity model which 
depends upon three Weiss fields &0j b\, 62 and involves 
three self-energy functions. When written in matrix form 
the impurity model must thus involve a closed algebra of 
three orthogonal symmetric matrices Mo, Mi, M2. Wc 
have not found a suitable closed algebra of 3 x 3 matrices; 
however a closed algebra of 4 x 4 matrices exists, with 



= 1, 



M 

M : = 

M, = 



1 
T2 



/0 1 IN 
10 10 
10 1 

Vi 1 0/ 



/0 1 0\ 

1 

10 

vO 1 0/ 



(27) 
(28) 

(29) 



This choice corresponds to a four-site real-space clus- 
ter with the same topology as that considered by Licht- 
enstein and Katsnelson in a pioneering study of su- 
perconductivity and antiferromagnetism in the Hubbard 
models 

The impurity model action becomes 



5(4) = Trm 



S 



(4) 



(30) 



leading to the self consistency conditions 
ISIylZ^ f 

l^— = J {dp)Mv)G(p,u) (31) 
with [z is the number of nearest neighbors.) 

0o = 1, (32) 



?i = e 



ip-a \ A 



t>2 = e 



i(p-a-\-p-b) 



,i(p-a+p-b) 



z(z-l) 



(33) 



(34) 



We finally discuss the interaction terms. In models 
with purely on-site interactions, the s\^ t may be fixed 
by analysis of the local limit, in which G a = for a ^ 
and Gq is known. For example, in the Hubbard model, 
Go has poles at u> = ±U/2. The absence of any intersite 
correlations ensures that any off-diagonal terms in S^j = 
and the identity of the n sites implies the interaction 
is just 



S 



(n) 



u 



j=l,...,n 



(35) 



When treated perturbatively to order U 2 , the two and 
four impurity models reproduce exactly the appropriate 
Fourier coefficients of the exact (perturbative) lattice self- 
energy of the Hubbard model. To gain some idea of ef- 
fects beyond perturbation theory we show in panels (a,b) 
Fig. 2] results for the real and imaginary parts of the self- 
energy obtained from the 'IPT' approximation, 1 in which 
the impurity model and self consistency condition are 
solved by writing the second order perturbation theory 
expression for the self-energy, but using exact (impurity 
model) Green functions. We observe that in the IPT ap- 
proximation to the Hubbard model, non local effects are 
very small. 



IV. CAUSALITY 

A. General considerations and a simple example 

A key difficulty in 'cluster' extensions of the coherent 
potential approximation or of dynamical mean field the- 
ory has been causality. Any physical theory must be 
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FIG. 1: (a) and (b): Real and imaginary part of self-energy 
obtained from 'iterated perturbation theory' approximation 
to half-filled Hubbard model (U = 4, T = 0) using one and 
four-site clusters. Dash-dotted curves are for one-site cluster. 
Solid, dashed, and dotted curves are So, £i, and £2 in four- 
site clusters, respectively. £1 and £2 are multiplied by 4. 
(c) and (d): Self energies in four-site clusters obtained from 
simple filtering procedure of section IV, with /1 = 0.85 and 

h = fl 



causal, which implies in particular that /mS(|),w) < 
0. Many cluster schemes however generate functions 
^approx(p, with the unfortunate property that for 
some range of p, u>, ImYl approx (ji,u!) > 0. These viola- 
tions of causality have apparently never been clearly un- 
derstood or cured, but are generally viewed a deficiencies 
of the 'cluster model' used to calculate the self-energy. 

The formal development of the previous section sug- 
gests that the causality violations may be thought of as 
an example of the 'ringing' phenomenon familiar from 
Fourier analysis. Any straightforward real-space cluster 
scheme corresponds to using some set of orthogonal func- 
tions to expand the momentum dependence of the exact 
lattice self-energy £(f>, uj) in the sense of Eq. ((TTJl . The 
first term in an orthogonal function expansion is <j>Q = 1 
which is everywhere positive, so the local approximation 
is guaranteed to be causal, but all of the other orthogonal 
functions change sign over the Brillouin zone, so an ex- 
pansion which is truncated at low order is not guaranteed 
to be positive everywhere in the zone. In particular, if at 
fixed u> the function S(p, uj) has a strong, narrow peak at 
some momentum p, then truncating an orthogonal func- 
tion expansion at a low order will produce a self-energy 
whose imaginary part changes sign. As noted above, in 
the "DCA" approach^ this problem is avoided, at the 
expense of introducing discontinuities in the momentum 
space representations of G and E, by choosing the or- 
thogonal functions to be a tiling of the Brillouin zone by 
rectangular filters. However, the general argument, that 
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FIG. 2: Momentum dependence of the second order perturba- 
tion theory approximation to the imaginary part of the self- 
energy of the half-filled two dimensional Hubbard model at 
T = 0. Solid, dashed, dotted and dash-dotted curves are for 
ui — 10, 2, 1 and —10, respectively. 

acausal behavior found in most schemes arises from 'ring- 
ing' associated with narrow peaks in the momentum de- 
pendent self-energy, suggests also that the phenomenon 
is not of very great significance. 

To illustrate the issue we show in Fig. |5]the momentum 
dependence of the second order perturbation theory ap- 
proximation to the imaginary part of the self-energy of 
the two dimensional square-lattice Hubbard model, for 
several different frequencies. At frequencies within the 
electronic band the momentum dependence of the self- 
energy is relatively weak, but at frequencies well above 
the upper band edge or well below the lower band edge, 
the imaginary part of the self-energy becomes sharply 
peaked in momentum space. This phenomenon has a 
simple kinematic origin. At this order of perturbation 
theory, the imaginary part of the self-energy at u) > 
corresponds to a decay of a particle into two particles 
and a hole. At u> = 10, energy conservation means the 
allowed final states correspond to two particles near the 
top of the band [momenta near (tt, tt)] and a hole near 
the bottom of the band [momenta near (0,0)]. Momen- 
tum conservation then restricts the initial momentum to 
be near (0, 0). 

We now consider expanding the self-energy in Fourier 
harmonics. The symmetry of the hypercubic lattice im- 
plies that 

E(p, w) = Eo(w) + 2dj p Ex(uj) +2d{d- 1) 7^ 2) £ 2 M + ■•• 

(36) 

with, for a hypercubic lattice of unit lattice constant in 
d dimensions, 

7 P = 2 cos (^)' ( 3? ) 

a— 1. . .d 

lp ] = J7jz 1) H cos(p Q )cos(p 6 ). (38) 

^ ' a=l...d 

The dashed and dotted curves in Fig. |3|show the re- 
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FIG. 3: Comparison of momentum dependence of exact sec- 
ond order perturbation theory approximation to self-energy 
of half-filled two dimensional Hubbard model (solid line) and 
the approximation using Eq. i'Mil at T — 0; dash and dotted 
curves are the approximation using the first and second, and 
first three terms in Eq. I|36^ . 



FIG. 4: Main panel: momentum dependence of Im E(w) at 
u = 9, 11 calculated for U = 1,2,4,6 by use of the 'IPT' 
approximation on a four-site cluster. The increase of U is 
seen to lead to a broadening of the peak, which for the large 
U values is seen (inset) to eliminate the 'acausal' behavior 
entirely. 



suit of approximating the exactly calculated £ by the 
first and second, or first three terms in the series given 
in Eq. An acausal behavior is observed at fre- 

quencies well above the upper band edge or well below 
the lower band edge, in accordance with the qualitative 
arguments presented above. At lo = 8, we observe the 
acausal behavior does not occur if only the first neighbor 
term is retained (dashed line) but does occur if both first 
and second neighbor terms are retained. This is because 
JmE 2 (w) changes sign between ui — 8 and 9 to fit the 
exact ImS(p, uj). [One can see this behavior even in Fig. 

HI(b,d).] 

As noted above the acausal behavior has (at least in 
this instance) a simple kinematic origin, which suggest 
that it may not occur in all circumstances. For example, 
at larger U, more complicated decay channels (e.g. one 
particle decaying into 3 particles and 2 holes) with fewer 
kinematic constraints may become important, leading to 
a broadening of peaks in the momentum dependent self 
energy and therefore to a causal self energy. An exam- 
ple of this phenomenon is shown in Fig 01 which displays 
the imaginary part of the self energy calculated high fre- 
quencies in the 'IPT' approximation at several different 
[/-values. The broadening of the peak and the destruc- 
tion of the acausal behavior suggested by the arguments 
above are clearly observed. We stress that the IPT is 
an uncontrolled approximation — however it is a compu- 
tationally tractable example illustrating a phenomenon 
which we suspect is of more general significance. 



B. Filtering 

We have argued that any cluster extension of the dy- 
namical mean field method amounts to a scheme for com- 
puting the coefficients in a low-order orthogonal func- 
tion expansion of the electronic self-energy. This sug- 



gests that the reported causality difficulties are generic 
and inescapable: no truncation of an orthogonal func- 
tion expansion is guaranteed to lead to an everywhere 
non-negative approximant for an arbitrary test function. 
We emphasize, however, that the impurity model itself 
is causal; it is only the resulting lattice self-energy which 
may have acausal features. 

The discussion of the previous subsection suggests that 
the acausal behavior is 'ringing' is associated with a rel- 
atively large peak near a particular momentum. It is 
possible that such a peak could be physically important, 
occurring for example at a fermi surface 'hot spot' in a 
system close to a quantum critical point. In such a case, 
long ranged interactions are evidently physically crucial, 
and modeling them with a local model is simply inap- 
propriate. However, as seen in the explicit calculations 
present above, peaks in JmS which are both sharp and 
large relative to typical values of £ are more often asso- 
ciated with band edges and extremal frequencies. These 
regions of u and p are not particularly important for en- 
ergetics, suggesting that the acausality is simply a minor 
technical annoyance. In what follows we consider meth- 
ods of removing it. 

We first note that in particular temperature regimes of 
particular problems (for example, the double exchange 
model on the simple cubic lattice and at not too low 
temperatures^) the ringing phenomenon might not oc- 
cur. We also note that a clever choice of expansion func- 
tions may mitigate the severity of the problem. For ex- 
ample, we saw in the Hubbard model the difficulties arise 
from narrow peaks centered near the band edges (p « 
and p ~ (tt, 7T, ..)). An expansion based on the functions 
1,7 P and a function orthogonal to both 1 and j p but 
strongly peaked near and (7r,7r, ..) might have a wider 
range of applicability than the straightforward harmonic 
expansion. 

A more general approach is to filter the self-energy, for 
example by convolving it with a function to smooth out 
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FIG. 5: Comparison of filtered, unfiltered and exact re- 
sults for second order perturbation approximation to self- 
energy of the half-filled two dimensional Hubbard model at 
u> = 10. Solid curve is the exact result. Dashed and dot- 
ted curves are obtained from four-impurity model with filter- 
ing coefficients /i = 1 and 0.85, respectively, and /2 = ft- 
Numerical calculation gives the following values for ImTn: 
(JmEo,7mEi,/mE 2 ) = (-0.00418,-0.00234,-0.00114) for 
/i = 1 and (-0.00418, -0.00144, -0.00033) for ,fi = 0.85. 
Calculations are done at T = 0. 



any sharp peaks and then approximating the smoothed 
function by a low order harmonic expansion. Let us make 
this approach more precise, writing Eq. (| 1 11) as 



J approx 



E 

i— 0..n 



fi<j>i(j>)Y,i{u), 



(39) 



where /o = 1 and 1 > /j>o y are the Fourier compo- 
nents of the smoothing function. Carrying through the 
development of the previous section leads to 



SCI 



approx 



<5Ej(w) 



= fi / (d P )G p (u J )Mp) 



(40) 



so that the quantum impurity model acts to reproduce 
the smoothed Green's functions of the theory. (We note 
that sum rules typically constrain the large u> behavior 
of the local Green function, so that one must choose 
/o = 1 in order to have a consistent representation at 
least within the simple impurity models of which we have 
studied). Fig. |5] shows the results obtained using differ- 
ent choices of simple filtering coefficients in the low order 
self energies discussed in the previous section. 

The straightforward filtering approach gives up some 
fraction of the intersite correlations in order to guarantee 
a causal theory. A more sophisticated possibility would 
be a frequency dependent filtering. We saw from the 
low order perturbation calculation on the Hubbard model 
that the self-energy was only very strongly peaked in mo- 
mentum space for very high or very low energy states, 
where decay kinematics constrained all states involved to 
be near the band edges. These states are not very impor- 
tant to energetics. One may therefore filter out only these 
by setting an arbitrary frequency scale uif and A and 
writinge.g. /i = 1— f/[exp{(— |w|+w/)/A}+l]. We note 




FIG. 6: Real (a) and imaginary (b) part of self-energy, Ei 
and E2 , obtained from 'iterated perturbation theory' approx- 
imation to half-filled Hubbard model (U = 2, T = 0) using 
four-site cluster with frequency-dependent filtering. Solid and 
dashed curves are the results with parameters (/, w/,A) = 
(0.15,6,1) and (0.15, 4, 1), respectively. Dotted curves are 
the results without filtering. 



that this filtering may be performed ex-post- facto: one 
may solve the impurity model, determine the band-edge 
regions where filtering is required, and then re-solve the 
impurity model with filtering only in these regions. Fig. 
shows the results obtained using the above frequency- 
dependent filtering with different choices of parameters. 
Here, we take fa = , and four- impurity model is self- 
consistently solved. It is clearly shown that £1,2 are re- 
duced at higher frequency region where acausal behavior 
has been observed in the second order perturbation the- 
ory approximation, while there is not much effect in the 
low frequency region. 

A still more sophisticated approach would be to mix 
different harmonic components, so that low order terms 
were modified only if higher order terms were impor- 
tant (i.e. if filtering were really needed). This approach 
is most helpful if a-priori knowledge of the locations of 
sharp momentum-space structures is available. Indeed, 
the calculations presented above of the behavior of the 
Hubbard model strongly suggest that the singularities are 
associated with the top and bottom of the band, i.e. with 
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states around p = (0, 0, ..) or p = (it, 7T, ..). To define this 
transformation we begin from the harmonic expansion, 
Eq. iJHJ, and then introduce a transformation which 
mixes the harmonic coefficients, so that 

T, approx (p,uj) = ^2 0t(p)/y£j(w)> ( 41 ) 

i,j=0..n 

and the self consistency equation becomes 

Gi(u) = J2 fa J ( d P) Gp(«)^(p). (42) 

3 

We have however not yet explored these more compli- 
cated filtering procedures 

V. CONCLUSION 

In this paper we have argued that one should take a 
more abstract view of cluster extensions of the dynami- 
cal mean field method; regarding them as algorithms for 
computing frequency dependent coefficients in an orthog- 
onal function expansion of the electron self-energy. This 
approach renders moot the debates about 'correct' choice 
of cluster and embedding, and clarifies the meaning of the 
'causality violation' encountered in real-space-cluster ex- 
tensions of the dynamical mean field method. We have 
suggested that the causality violation is in most cases a 
minor technical problem which can be cured if needed 
by any of a variety of 'filtering' procedures. 

Several questions arise, for which further research 
would be desirable. One is the question of the cor- 
rect choice of interaction terms in the "fictive impurity 
model". This issue has been discussed in the context 
of on-site interactions, especially in connection with the 
"DCA", 14 and, for longer ranged interactions in the con- 
text of the E-DMFT approachiSiiSi but deserves more at- 
tention in the context of the more conventional orthog- 



onal function expansions. A second key issue is the ori- 
gin of the 'causality violations'. We have shown that 
in simple perturbative models (chosen because exact ex- 
pressions for the momentum dependent self energy are 
available) the violations are an example of the famil- 
iar 'ringing' phenomenon, and moreover occur mainly in 
band edge regions of little kinematical importance. The 
issue however deserves further exploration in less triv- 
ial contexts. A third open problem is the question of 
which approximants to the self-energy are representablc 
by impurity models. We have presented arguments (sub- 
stantiated by low-order perturbative calculations in dif- 
ferent limits) indicating that straightforward orthogonal 
function expansions, with and without filtering, are rep- 
resentable. However, one might imagine that more com- 
plicated approximate representations of the self-energy 
might be advantageous in some problems. We do not 
at present have a general framework for determining the 
circumstances under which a general approximant to S 
is representable in terms of an impurity model. The issue 
of different choices of filtering function seems also likely 
to benefit from further research. In particular, it seems 
likely that one can more efficiently exploit the observa- 
tion that the acausality is associated mainly with the 
p = and p = (tt, 7T, ...). Finally, applying the method 
to a wider range of models, to explore which choices lead 
to good representations of the physics of interest, is an 
urgent task. 
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